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Abstract 

We compute the static energy of QCD at short distances at next-to-next-to-next-to leading- 
logarithmic accuracy in terms of the three-loop singlet potential. By comparing our results with 
lattice data we extract the value of the unknown piece of the three-loop singlet potential. 
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I. INTRODUCTION 



The energy between a static quark and a static antiquark separated at a distance r, which 
is referred to as the static energy, is a basic object to understand the dynamics of QCD 
When calculated at short distances in perturbation theory, the virtual emission of ultrasoft 
gluons (i.e. gluons with energy and momentum smaller than 1/r that can change the color 
state of the quark- antiquark pair from singlet to octet) produce infrared divergences. These, 
after resummation of certain diagrams, induce in the static energy logarithms, lnas(l/''") 
In an effective field theory framework 0, 4], the separation of scales in the problem is made 
explicit: the static energy becomes the sum of a matching coefficient, the static potential, 
which encodes all contributions from the scale 1/r, and of the contributions coming from 
ultrasoft gluons, which start at three loops. Logarithms are remnants of the cancellation 
between infrared divergences of the potential and ultraviolet divergences of the ultrasoft 

ri 

contributions 5]. These logarithms may be potentially large when r is very small, so that 
Q;s(l/r) lnQ;s(l/r) terms in the static energy may need to be resummed. The resummation 
of the ultrasoft leading logarithms (LLs)^ was performed in 6j. The anomalous dimension 
of the ultrasoft next-to- leading logarithm (NLL)^ was calculated in 3]^ and turned out to be 
quite large, showing that the resummation of the ultrasoft NLLs should also be addressed. 
In this paper, we will consider this resummation. 

n n n n 

It has been arg ued in [9|, M, lUl, UM that the proper consideration, and cancellation. 



of the renormalon singularities is crucial to obtain a good convergence of the perturbative 
series for the static potential in the short-distance region. The detailed analysis of the 

the potential will not be 
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121 ]. which essentially 



possible influence of ultrasoft effects in the renormalon structure o: 
presented in this paper. We will just follow the analysis of 
take advantage of the fact that a constant term may be added to the potential, a freedom 
that remains even if ultrasoft effects are taken into account. As we will see later, this 
seems to be enough to obtain a convergent perturbative series for the static potential in the 



^ Leading ultrasoft logarithms contribute to the static energy at order a^^" In" for n > 0, i.e. at next- 
to-next-to-leading logarithmic, N^LL, order. 

^ Next-to-leading ultrasoft logarithms contribute to the static energy at order a^^" In" for n > 0, i.e. at 
next-to-next-to-next-to-leading logarithmic, N^LL, order. 

^ The calculation relies on the next-to-leading order (NLO) calculation of the chromoelectric correlator done 
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short- distance region. 

The static potential is a basic ingredient of heavy-quarkonium physics |13]. In particular, 
its perturbative evaluation at higher orders is relevant to describe the top-quark pair pro- 
duction process near threshold. This process is expected to allow the extraction of the top 
quark mass to a high precision, and hence a remarkable effort is being made to calculate it 
at next-to-next-to-next-to- leading order (N^LO) ll5|, [ly, [iTj] . The past experience with 
the next-to-next-to-leading order (N^LO) results ^ indicates that both renormalon cancel- 
lation and the logarithmic resummation 19|, |20| are necessary for accurate determinations 



of the position of the pole and of the shape of the cross section respectively. Our results will 
be relevant for the N^LL calculation of this process. 

The paper is structured as follows. In the next section, we introduce residual mass 
terms in potential Non-Relativistic QCD (pNRQCD) and summarize the current status 
of the perturbative calculations for the static potential. In section llllt we present and 
solve the renormalization group equations for the static pNRQCD Lagrangian, at next-to- 
leading order, and briefly describe the renormalon subtracted scheme that we use. Section 
IIVI presents a comparison of our results for the static energy with lattice data, and the 
numerical extraction of the missing piece of the three-loop static potential. We conclude in 
section |Vl 

II. POTENTIALS AND RESIDUAL MASS TERMS IN PNRQCD 

The general form of the dimension 6 operators in the pNRQCD Lagrangian is 

CsS^S + CoO^^O", (1) 

where S is the singlet and 0° the octet fields. The coefficients Cg and Cq have dimension 1. 

Let us recall that, in order to define Heavy Quark Effective Theory (HQET) beyond 
perturbation theory, or even in perturbation theory when regularizations with an explicit 
scale (cut-off) are used, one needs to introduce a residual mass term 6mQ in the Lagrangian 



2l| 



(2) 



CuQET = {tDo - 6mQ) ^ + { — ] , 

with rriQ the heavy-quark mass and ip the heavy-quark field. We may associate to Sniq the 
size of the typical hadronic scale Aqcd- This residual mass term will be inherited by the 



pNRQCD Lagrangian. In the paper, we consider the weak-couphng regime of pNRQCD, 
defined by 

1 

(3) 



1 CKg . 

- > — > Aqcd, 



at leading order in the l/mg expansion; in this situation, the residual mass term is absorbed 
in the coefficients and Co above. Therefore, it is useful to split them in a part that is 
proportional to 1/r, which corresponds to the singlet and octet potential^, and a part that 
is proportional to Aqcd^ 



Cs ^ Vs+ As 
Co = Vo + Ao 



r 
r 



(4) 
(5) 



where 



, _ as(l/r) 



+0^2 s,o 
+ 



«s(l/r) 



+ 



16 vr^ 



47r 



«s(l/r) 

471 
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aPln^r/i+ ( - ^n'^ C\Pq{5 - 6 In 2) ) lnr/i + 04^,0 



with 



as(l/r) 
47r 



+ 



(6) 



ai 

a2,s 



31 10 



/4343 
VT62" 



|-8C(3))c.n,+ (^n, 



TT 



899 

sT 



+ ( y - 47I) /3o + (4ai/3o + 2/3i) , 



167r2^3 / 11^ 2 , 
--C. + -n,), 



ai + n/(-- + -ln2)+C^ 



149 22, ^ 4 2 

In2 + -7r2 

27 9 9 



(7) 



(8) 

(9) 
(10) 

(11) 



* The singlet potential is often referred to as the static potential, a terminology which we also adopt in 
the paper. Recall that it coincides with the static energy up to two loops but differs from it beyond that 
order. 
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/i is the ultrasoft factorization scale. The color factors are defined as Cp = {N^ — l)/{2Nc), 
Ca = Nc where Nc is the number of colors; ^/ is the number of (massless) flavors; 'Je is the 
Euler constant. The strong coupling constant ttg is in the MS scheme. The beta function is 
defined as 

7 9 °0 

«s/3(«s) = :nr^ = -:eEfc (12) 



dliifi 

where po = llC^/3 - 2n//3, A = 34^2/3 - 
of (32 and (3^ may be found, for instance, in 



n=0 



22 



OCArif/S — 2Cpnf, and exphcit expressions 
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The one-loop coefficient Si was calculated in [2J, |25|] , the two-loop singlet coefficient 02 



m 



29| and the two-loop octet coefficient 



the third-order correction was calculated in 



a^s,o has not been complete 



Q,|5|,Q^ 



00,0 m 



30|. The logarithmic piece of 



31| , whereas the non-logarithmic piece 



y calculated yet. The fermionic contributions of a^s has been 
presented very recently in [3^, where the computation of the nj independent piece is reported 
to be in progress. A Pade estimate of a^^s gives: 0.^1.,^ = — 487r'^ vi^'' , Vs'^\nf = 3) = —38.4, 



vPijif = 4) = -28.7, V^-'^rif = 5) = -20.5 |33|. The double logarithmic coefficient 



K3) 



may be obtained from and the logarithmic coefficient was obtained in [7|^. h.s,o 

stands for As^o(r, /i)^. 

At order in the multipole expansion, the dimension 6 operators of pNRQCD do not 
have an anomalous dimension and, therefore, the renormalization group equations for the 
coefficients A^ and Aq will have the same structure as in the HQET the renormalization 
group equation for the coefficient of the operator ip'^ip has. At next-to-leading order in the 
multipole expansion, the pNRQCD Lagrangian reads 



-pNRQCD 



light 



dh Tr <^ \id< 



CAT 



,^)]S + Ot [zDo-c„(r,^)]0 



+14 (r, ^)Tr {OV ■ S + SV ■ ^E 0} 
{OV • (7E O + QtOr • gY.) 



(13) 



where £iight is the part of the Lagrangian involving gluons and light quarks, which coincides 
with the QCD one, S = IJy/WlS, O = V2T"0", E is the chromoelectric field, Va and Vb 
are matching coefficients associated with the 0{r) operators of the pNRQCD Lagrangian 



^ Only the coefficient of tlie singlet potential was obtained there, it will be shown later in the paper that it 

coincides with the coefficient of the octet potential. 
® To simplify the notation, we will often suppress the dependence on r and just write As.o(/i). 
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and the dots stand for higher-order terms in the multipole expansion. Ultrasoft gluons cause 
transitions between singlet and octet fields and generate an ultrasoft anomalous dimension 
for the dimension 6 operators. In particular, this modifies the renormalization group (RG) 
equations for and A^. 



III. RENORMALIZATION GROUP 



The general structure o: 
case has been discussed in 



the renormalization group equations of pNRQCD in the static 
6|, where the complete N^LL order was calculated. Here, we will 



calculate the complete N^LL order. It has been proved in ^| that in order to perform the 
calculation one needs not to consider higher orders in the multipole expansion beyond those 
already contributing to the N^LL calculation. Hence, the structure of the RG equations 
remains the same as in [Q] , but the anomalous dimensions need to be calculated to one order 
more in the ultrasoft loops. The RG equations read 



d 

dn 



d 



3 „2 



'ys{as)VJ^ (c^ - c^) r 



d 



dfi 



(14) 



d 

djj, 
d 

dfi 

where the anomalous dimensions 7s(as), 7o(«s), 7a(«s) and 7_b(q;s) are needed at order a^. 
Strictly speaking the equations above hold for Cs,o = K,o + ^s,o, provided that Kj,o ^ ^s,o 
and one stays at linear order in Ag^o- If quadratic or cubic terms in A^ o are included, 
additional counterterms in the potential are needed to absorb the ultraviolet divergences of 
the ultrasoft calculation. 
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A. Anomalous dimensions 



We calculate the anomalous dimensions in a regularization scheme in which the gluons, 
light quarks and center of mass motion are taken in D dimensions but the potentials in 
the ultrasoft loops are kept in three dimensions. We renormalize using the MS scheme in 
(relative) coordinate space. 

In this scheme, the anomalous dimension 7^ was obtained in [3]. It is — xdZ^^^ /dag, 
where Z^^^ denotes the coefficient of the 1/e (1/e = 1/e — + ln47r, D = A — 2e) pole of 
the ultrasoft contribution to the static energy; we have 

7.(a.) = -|— fl+6-B). (15) 



3 TT V TT 

with 

5 ^ -5nf + CAjQn' + 47) 

108 ■ ^ ^ 

In a similar way, 70 can be obtained from the 1/e poles of the ultrasoft contribution to 
the self-energy of the octet field. It is important to recall that, although we may obtain it 
by matching gauge- dependent Green's functions, the self-energy of the octet field is a gauge 
invariant quantity in perturbation theory, for the same reason as the pole mass is. At the 
order we are interested in, it can be obtained from the following expression to be taken in 
the T 00 hmit fl 



2jVc J-T/2 J-T/2 



X (0(T/2, t):Jr ■ E«'(t)r ■ E"' {t')<P{t' , -r/2)^,'S)) . (17) 

T°-WciT^ stands for a and a insertion at the time T/2 and — T/2 respectively in 
the space sides of a rectangular Wilson loop, and 4>(t,t')^^ is a Wilson line in the adjoint 
representation. The Feynman diagrams involved in the evaluation of this quantity in a 
covariant (or Coulomb) gauge do not coincide with the ones needed for 7^ (compare Fig. [T] 
with Fig. 4 of 0]), and would require extra calculations. Fortunately there is an argument 
which makes the explicit calculation unnecessary. If we take the Aq = gauge, the number 
of diagrams to be evaluated collapses to a few (the octet field does not emit gluons anymore), 
which, in addition, are the same for 7^ and 70. In this gauge, both anomalous dimensions 



are related by trivial color factors. Since both anomalous dimensions are gauge invariant, it 
turns out that we can read 70 from the known result for 7^: 

7oK) = (18) 

At order as, this is confirmed by the explicit calculation of {4]. 

Let us turn now to the evaluation of 7^ and 75. We use the fact that the /i dependence 
of Va and Vb can be also obtained from the infrared logarithms in the matching calculation 
between HQET and pNRQCD. There is only one diagram which is infrared divergent in the 
matching calculation: it is displayed in Fig. [2l However, the divergence turns out to be 
linear and produces no logarithms. Then these anomalous dimensions remain zero also at 
next-to-leading order: 

7^(as) = 7b(«s) = 0. (19) 



B. Solution of the renormalization group equations 

The renormalization group equations for Vg and Vq can be read from 0141) using that 
Vs^o 3> As,o and neglecting the latter. They are given by 

3 „2 





2a,C^ 

3 TT 


d ^ 
dfi " 


las 1 


d 

djj, 




dfi 





-Vb = 
dfi 






(l + Q'^B)vl{Vo-Vsfr 



\3 „2 



(20) 



The solutions of these equations, at the order we are interested in, are Va{ij) = Va(1/t) 
Vb{^) = VB{llr) = 1 (from VaA^I^) = 1 + 0{al), see fl), and 



K(/i) = Vs{l/r) + -CFr^[Vo{llr)-Vs{llr)f 



(a) (b) 




(k) (1) 



FIG. 1: Ultrasoft contributions to the self-energy of the octet field in a covariant gauge, at the 
order of interest. Symmetric graphs are understood for (c)-(i). 
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FIG. 2: Infrared divergent diagram in the matching calculation of Va and Vb- 

5l7^ ^ * -"'''^'■")' '''' 

VM = K,(l/r) - gi-r^ |\/„(l/r) - V.(l/r)|' 

where 

1 / ft 12B\ 

The renormahzation group equations for and Aq can also be obtained from f|T4l) by 
expanding c^^o about K,,o and keeping the terms hnear in Ag o- They are given by 

' f^^As = (1 + 6^5) Vy [K{l/r) - K(l/r)]' (A, - A,) 

an 71 \ TT / 

Ct/i 

where we have aheady approximated and by Vo{l/r) and Vs(l/r) (the /i de- 

pendence of Vo{fi) and enters at N^LO, which is beyond the accuracy of ( |24l) ). The 

solutions of the renormahzation group equations read, 

A,(,^) = A',A + 2CF(A'„-A'.)Ar^|V;(l/r)-K(l/r)l' 
A„(,,) = A'„A--i(A'„-Af.)Ar'|K,(l/r)-K(l/r)]' 

where Ns, No are two arbitrary scale- invariant dimensionless constants and A is an arbitrary 
scale-invariant quantity of dimension one. 
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The integration constants A^^ and No are fixed by the initial conditions, As(l/r) and 
Ao(l/r), of the solutions of the RG equations. In turn, the initial conditions are fixed by 
matching a suitable Green's function in QCD with the corresponding one in pNRQCD. 
Note that if at the matching scale As{l/r) = Ao(l/r) = 26mQ, as it happens in MS-type 
schemes, then As(yu) = Ao(yu) = 26mQ for any fi. We will see, in the following sections, the 
convenience to use an RS (renormalon subtracted) scheme. In the RS scheme As(l/r) and 
Ao(l/r) are different constants (which also differ from 25mQ) and hence As(/i) and Ao(/i) 
evolve in a non-trivial way according to the RG equations above. 



C. The Renormalon Subtracted scheme 



The discussion in the previous sections is independent of the renormalization scheme 
used for the matching calculation between HQET and static pNRQCD at the scale 1/r. 
The outcome of the matching calculation only enters through the initial conditions of the 
RG equations. It is well known that the singlet potential Vg calculated in the MS scheme 
displays a bad behavior as a series in as(l/''") even at small values of r. This bad behavior may 
be ascribed to renormalon singularities that lie very close to the origin of the Borel plane. In 



order to treat the renormalon singularity, we shall follow the procedure described in Ref . [3J] , 
the so-called renormalon subtracted (RS) scheme. Under RS scheme we understand a class 
of subtraction schemes that subtract from the perturbative series of Vs^o in the MS scheme 
the non-integrable piece at u = 1/2 in the Borel transform of the potential, expanded about 
u = and integrated over u. The whole non-analytic piece (non-integrable and integrable) 
at M = 1/2, expanded about m = and integrated over u, reads (at the scale p) 

n=l ^ ^ fc=0 ^ ' 

The coefficients are given in terms of the coefficients of the beta function. Since the 
beta function is known up to four loops only, all d^ for A; > 3 are unknown; the known terms 
are 

c/o = 1, 

P>1 - /32/3o 



-2/3o^/33 + 4/53/3i/52 + iil m - 2Pf) - 2(5^^2 + ii\ ^^^^ 



32(6 - \)hi5l 
11 



(29) 



with 

Hence, it is practically unfeasible to subtract Eq. fl271) . This is not a real problem because 
only subtracting the k = term, which corresponds to the non-integrable piece in the Borel 
transform, is necessary in order to obtain a series that is Borel summable. In the following, 
we will subtract all the known terms in Eq. fl27|) . i.e. up to A; = 2, as was done in the 



original proposal of the RS [341] • 



Furthermore, since the potential is given as an expansion of as(l/r), in order to achieve 
a successful renormalon subtraction at every order in ag, it is important to expand as(p) in 
terms of as{l/r) (or viceversa). We chose to expand as{l/r) in terms of as{p) in Eq. ([6]) 
instead of doing the reverse in Eq. ( |2711 . because the uncertainty in the normalization con- 
stants Rs^o of the renormalon singularities is then largely absorbed in the arbitrary additive 
constant needed to compare with lattice data (as it will be described in the next section). 
This expansion generates Inrp terms, which will be kept from becoming large by choosing 
p = l/{r) = 3.25/ro, (r) being the central value of the range where we compare with lattice 
data and ro being the reference scale used in the lattice computation (see the next section). 
At the order we are working, we only need to keep terms up to order (Xg{p). 

IV. COMPARISON WITH LATTICE RESULTS 



In this section we will compare our results with the [uf = 0) lattice data of Ref. 35] . 
This will allow us to extract a value for the three- loop coefficient a^g. 

A. Setting the scales and parameters 

We choose, as anticipated in the previous section, p = l/(r) = 3.25/ro (the reference 



scale ro, used in the lattice computation, has a value of about 0.5 fm, see [35|| for more 
details; we will present all our results in units of ro). The remaining scales and parameters 
entering in the expressions are chosen as follows. The number of light flavors rif is set to 
zero. The ultrasoft scale /x is set to /i = 2/ro. The normalizations of the u = 1/2 renormalon 
singularities for the singlet and octet potentials, Rg o, are determined using the procedure 
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described in [36!]; one obtains: 



Rs 

Rn 



-1.333 + 0.499 -0.338 
0.167- 0.0624 + 0.00972 



-1.172, 
0.114. 



(30) 



as at the relevant scales is determined according to [37|, which uses Aj^tq = 0.602(48). We 
will use the running of according to the order we are working at (for instance, for the 
one-loop curve we use the two-loop running for the order as term and the one-loop running 
for the order a^ term, and so on). 

The three-loop coefficient a^ g, which enters in our N^LL results, is unknown. We can use 
the Pade estimate of |33| to get an idea of its expected size (that is, use Cq^^=o = 313 from 

~ (3) 

Table 1 of that paper, which corresponds to 03 = 114633; the relation between cq, Vs and 
03 s is given in Eq. (29d) of 33|], for easier reference we reproduce it in the appendix). To 
estimate the uncertainty that we should associate to that Pade value we can make use of 
the results of 32 1. In Ref. j32j], all the fermionic contributions to the three-loop coefficient 
Co are calculated, therefore the difference 



Co,nf ~ Co ,nf — l 



f- 



(31) 



is known (co,n/ stands for the coefficient Cq calculated with Uf flavors). Since Ref. ^] 
presents the Pade estimated values of Cq for nj = 0, ... ,6, we can check if those results 
satisfy the known values for (!3T|) or not. This comparison is presented in Table In 
addition, we can also obtain co for Uf = 0, which is the coefficient we need here, from the 
Pade estimated value of co for ri/ = 6 and (|3Tll :^ 



c, 



from n 



0,n/=0 



239, 



(32) 



to be compared with the value co,nj=o = 313 presented in [331]. We will take this as an 
indication that one should assign an uncertainty of around 30 — 40% to the Pade estimate. 
As we will see later, though, the lattice data is precise enough to be used to obtain an 
independent extraction of the value of Cq. 



^ In 32], the coefficient of the d' 



'ahcd ^ahcd 



color structure is given numerically, but the limited numerical 



precision is not yet affecting the numbers presented in Table ID 

We chose Uf — 6, because it is for Uf — 6 that the Pade approximation comes closest to the three- loop 
RG accessible coefficient Ci, in the notation of [sS^ that we reproduce in Eqs. (jA2[) and (IA4I) . 
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Exact result [32] 


Fade estimate [33] 


Uf 

J 


= 6 


-21.39 


-29.6 




= 5 


-26.56 


-37.4 


^/ 


= 4 


-31.86 


-44.5 


nj 


= 3 


-37.28 


-51 




= 2 


-42.84 


-57 


Uf 


= 1 


-48.52 


-63 



TABLE I: Values of the difference co^nf — co.n^-i for the exact result and the Fade estimates. 

B. The static potential 

Using the choices of scales and parameters described in the previous section, we obtain 
the singlet potential in the RS scheme represented in Fig. [3] (to alleviate the notation we 
do not explicitly indicate the dependences on /i and p of the different functions in the labels 
of the plots). In all the plots in this section, the dotted blue curve will be at tree level, 
the dot-dashed magenta curve will be at one loop, the dashed brown curve will be at two 
loop plus N^LL resummation and the long-dashed green curve will be at three loop (with 
the Fade estimated value for cq, i.e. cq = 313) plus N^LL resummation. As expected and in 
sharp contrast to what would happen in an on-shell scheme 13] , we see that when we use a 
threshold scheme that cancels the leading renormalon, like the RS scheme, the perturbative 
series for the potential exhibits a convergent behavior. 

C. The static energy 

In order to perform a comparison with lattice data, we have to plot the static energy Eq 
as a function of r: 

Eo{r) = Vs + As + 5us = K (r, p, p) + Ki{p) + K2{p)f{r, p, p) + (5us (r, p) , (33) 

where is given by Eq. fl25l) and 5us contains the contributions from ultrasoft gluons. 
Ki and K2 have to be understood in the RS scheme, which is where the p dependence comes 
from (at the order we are working, / will not depend on p, which will be dropped from it 
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FIG. 3: Static potential roV^(r), in the RS scheme, as a function of r/vo- The dotted blue curve is 
at tree level, the dot-dashed magenta curve is at one loop, the dashed brown curve is at two loop 
plus leading ultrasoft logarithmic resummation and the long-dashed green curve is at three loop 
(Fade estimate) plus next-to-leading ultrasoft logarithmic resummation. The solid black curve is 
also at three-loop plus next-to-leading ultrasoft logarithmic resummation but using cq = 250.7, see 
section IIVDI 

in the following). Also, at the order we are working, the renormalized expression of ^us is 
only needed at leading order^ (and can be read from equation (14) of [t^): 



6 



us 



C.fi^ai(l/r) 

24 r 71 



-2 In 



«s(l/r)iVc 
2r fi 



- - 21n2 
3 



(34) 



Ki and K2 are the constants A^'^A and {No — Ns)A in Eq. while f{r,fj,) is 
/(r, /i) = 2CFr' [Voil/r) - V,{l/r)f In + r/o [«s(/i) - «s(l/r)] 



(35) 



We are considering the weak-coupling regime in the static limit, defined by the hierarchy 
of scales given in ([3]). To have a definite way to organize the terms in Eq. (|33l) we will also 
use 



a: 



A ~ A'.A ~ NoA ~ Aqcd ^, 



r 



(36) 



^ Note that large logarithms have been resummed in Vs, so that the counting of (5us is a fixed order one if 
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^ ^ ^ , , ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ J 

0.0 0.1 0.2 0.3 0.4 

r/ ro 

FIG. 4: Illustration of the hierarchy of scales relevant for the static energy. The scales are: 1/r 
(dotted red), (Vq — Vs) (dot-dashed green), as{Vo — Vs) (dashed black) and = 0.602/ro (solid 
blue). {Vo — Vs) is taken at tree level, {Vo — Vg) = Ncas/{2r), and the one-loop running of 
as = as(l/r) is used. 

which is compatible with Eq. ([3]). Figure H] shows that the scale hierarchy of Eq. ([3]) as well 
as the counting above hold for t/tq < 0.5. 



The lattice data of 351] is presented as the difference between the static energy at distance 
r and the static energy at a reference scale r^. Therefore, what we actually need to plot, in 
order to compare with the lattice data, is 



Eo{r) - Eo(rmi„) + ^^oM^min) =Vs + Ki + K^f + 5us, (37) 

where rmin is the shortest distance at which lattice data is available, to-Eq'^**' (rmin) = —1.676 
is given in Table 2 of 35], and Ki is a suitable constant, which can be obtained by imposing 
the equation to be true at r = rmin. 

From the counting described above, we see that the last two terms in the rightmost part 
of Eq. fl?r|) only start contributing at the three- loop level. We also see that, at three- loop 
level, we only need the function /(r) in Eq. fl35l) at leading order, i.e. 

/(..)^2C.(^ + C.a.(lA))^|;..>^. (38) 
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Therefore, the static energy (that we will use to compare with lattice data) at N^LL order 
is given by Eq. fl37|) with Vs given by fl2T]) (understood in the RS scheme), / given by fl38|) 
and 6\js given by 0341) . The constant K2 will be fixed by a fit to the lattice data below 
r/ro = 0.5. In all the plots, we will always display our results until r/rg = 0.5, which is the 
region where we expect perturbation theory and our hierarchy of scales to be reliable. 

1. Analysis of the uncertainties in the static energy 

As we have already mentioned, our N^LL results depend on the unknown three-loop coef- 
ficient Co (in addition to the constant K2). The static energy also suffers from uncertainties 
due to the parameter, used to determine as, and from uncertainties due to the neglect- 
ing of Os/r and higher-order terms. Since we want to use the lattice data to extract cq, we 
need to make sure that the static energy is more sensitive to variations of cq (and K2, which 
will be also fitted to the data) than to the errors due to Aj^ and the higher-order terms. 

According to the discussion in section IIV Al we will let cq vary by 35% around the Pade 
value Co = 313. will be varied from -5 to 5. Aj^ will be varied according to the range 
quoted in Ref. (37|, A^ro = 0.602(48). Finally we assess the impact of neglecting the 
four-loop, ctg/r, terms by simply adding the term ±10 Ci^a^/r to the three-loop curve (±10 
is intended to be a rough estimate of the size of the four- loop coefficient in the RS scheme) . 
All those variations are shown in figure [5] as the green bands. We have used the Pade value 
Co = 313 for the plots where we vary K2, Aj^ and the higher-order terms. K2 = is taken 
for the plots where we vary cq, A^^ and the higher-order terms. We can clearly see that the 
variations due to co and K2 produce larger bands than those due to the uncertainties in Ay^ 
and higher-order terms. Therefore, it will make sense to fit co and K2 to the lattice data. 

D. Lattice comparison and extraction of cq 

We are now ready to compare with lattice data. As we have seen, the N^LL expression 
for the static energy depends on two unknown parameters, K2 and co, which we will obtain 
from a fit to the lattice data points below r/ro = 0.5. The result of this two parameters fit 
(Aj^jgTo = 0.602 and higher-order terms are set to zero) gives the values 

K2 = -1.0465 , Co = 250.7, (39) 
17 
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FIG. 5: Impact of the variation of cq, K2, and effect of higher-order terms in the N^LL 

expression for the static energy respectively, represented as the green bands, cq has been varied 
by 35% around the Fade value cq = 313. K2 has been varied from -5 to 5. has been varied 
according to Aj^tq = 0.602(48) [37] and the effect of higher-order terms has been estimated by 
adding the term ^IQCpa^/r to the N'^LL curve, cq = 313 and K2 = have been used as central 
values. The red points are the lattice data of [s^. 

the x^/d.o.f. of the fit is 0.07 (we have four degrees of freedom). Those are the best fit 
values, but to obtain an allowed range of values for cq, according to the lattice data, we 
will analyze the of the N^LL curve for different values of the cq and K2 parameters 
(with Af^jgTo = 0.602 and higher-order terms set to zero). We choose the allowed range 
of values for cq by identifying the region of the CQ-K2 parameter space where the reduced 
of the N'^LL curve is better than that of the N^LL curve, and, at the same time, K2 
retains a reasonable power counting value. The reduced of the N^LL curve is 3383 (in 
this case we have six degrees of freedom). We will allow values of \K2\ up to \K2\ = 2, to be 
conservatively consistent with our power counting. Figure |6] presents a scatter plot which 
shows the values of the ratio {x^ /d.o.f.)NZLL/{x^ /d.o.f.)^^^ with different values of the K2 
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FIG. 6: ix^ /d.o.f.)]\;3iL/{x^ /d.o.f.)]\[2ii for different values of the K2 and cq parameters. The 
color of each point in the cq — K2 plane represents the value of the ratio, lighter colors correspond 
to higher values and darker colors to lower values. We show values of the ratio up to 1. 

and Co parameters. In that plot, lighter (darker) colors correspond to higher (lower) values 
of the ratio. From that we obtain the range (215,350) for Cq, as can be read from the figure. 
To make the values of easier to visualize, we also present, in figure [3, separate scans over 
values of cq and K2, i.e. we scan over cq or K2, fit the other parameter for each point, and 
present the ratio of x^s as a function of Cq or K2. 

In figure [HI we show the N^LL curve, with the best fit values for K2 and cq given in (l39ll . 
together with the lattice data. We also show in the plot the tree-level, one-loop and N^LL 
curves for the static energy. Several comments are in order. First, we can see that the 
agreement with lattice is improved when we go from tree level, to one loop, to N^LL. We 
also note that the N^LL curve describes very well the data. The fit of the N^LL curve to the 
lattice data was not constrained to give a value of K2 compatible with the counting, but this 
turns out to be the case (recall that our power counting requires 17^21 ~ «s/''^ ~ A ~ 0.6, 
in units of rg). This gives us much confidence in the consistency of our analysis. Note 
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FIG. 7: {x^ / d.o.f {x^ / d.o.f .)j^2i^]^ for different values of cq (left) or K2 (right), the other 

parameter is fitted. 

also that the best fit value for cq is smaller than the Pade estimate of [33], but in better 
agreement with (132|) . We would like to emphasize that the values of Cq in Eqs. (!32|) and (!39|) 
are obtained by completely independent procedures. Finally, let us also note that, when we 
use the value of Cq in Eq. (139|) . the convergence of the perturbative series for the potential 
seems to be slightly improved, with respect to using the Fade value co=313 (see the solid 
black curve in figure [3l) . 

In figure [9] we present the bands obtained by varying cq and K2 according to the ranges 
described above, and the bands induced by the variations in and higher-order terms. 
All the bands are obtained by keeping the rest of the parameters at their central or best fit 
values. 

Let us comment at this point on the scheme, factorization scale and implementation of the 
RS scheme dependences of the extracted value of Co (or equivalently a3,s). Concerning the 
scheme dependence, recall that the scheme of Ref. has been used to factorize the ultrasoft 
contributions, which differs from standard MS scheme (we note that ultrasoft contributions 
were ignored in the Fade estimate of Ref. 33|). However, if we add to 03 ^ the corresponding 



non- logarithmic parts coming from (5us; then we obtain the non- logarithmic piece of the static 
energy at three loops (which was denoted as 03 in j3]). The three- loop coefficient of the 
static energy 03 is independent of the scheme used to factorize the ultrasoft contribution, as 
opposed to 03 s. At the practical level, the non- logarithmic part of the ultrasoft contribution 
just shifts the values of Cq that we quote above by —1%. Concerning the factorization scale 
and implementation of the RS scheme dependences, we have checked that varying Rg by 
30%, Ro by 10% (according to the corresponding last known terms of the series in ( l30l) ). the 
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FIG. 8: Comparison of the singlet static energy with lattice data^ 
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ilatt. 



We plot 



35( 1 (red points). 



(^min)) a function of r/ro and the lattice data of 
The dotted blue curve is at tree level, the dot-dashed magenta curve is at one loop, the dashed 
brown curve is at two-loop plus leading ultrasoft logarithmic resummation and the solid black curve 
is at three-loop plus next-to-leading ultrasoft logarithmic resummation, using the best fit values of 
Eq. (I39D. 



scales p and p by 10% and implementing the RS subtraction with just the do term in fl27|) 
shifts the range for cq by a maximum of 7%. 

We conclude by noting that both the best fit value from the lattice comparison and Eq. 
f l32|) indicate a value of Cq lower than the Pade estimate Cq = 313. The computation of 
this three-loop coefficient is reported to be in progress 321]. For the sake of comparison, the 
static energy at N^LO is given by 



En(r) 



+a2 



Cpa,{l/r) \ . a,(l/r) 

i + 0-1 



«s(l/r) 



47r 



+ 



47r 

l^Cl In ^-^"-(^/^^ 



+ as 



as(l/r) 



47r 



(40) 



with a2 = and as in the range (1.08, 1.17) x 10^. 
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FIG. 9: Impact of the variation of cq, K2, and effect of higher-order terms in the static energy, 
respectively, represented as the gray bands, cq has been varied in the interval (215, 350), K2 in 



the interval (-2, 2), Aj-jg has been varied according to Aj-jg-ro = 0.602(48) 



37|] and the effect of 



higher-order terms has been estimated by adding the term ±10 Cpa^/r to the N^LL curve. All the 



bands are obtained by fixing the other parameters at their centra 
are the same as in figure [HI The red points are the lattice data of 



or best fit values. The curves 



V. CONCLUSIONS 



We have calculated the QCD static energy at short distances at N'^LL accuracy, in terms 
of the three-loop singlet potential, whose coefficient a^ g for Uf = is the only missing 
ingredient in our calculation. It is remarkable that such a higher-order calculation can be 
carried out analytically, which shows, once more, what invaluable tools effective field theories 
provide for higher-order calculations. The static energy at this order turns out to depend 
on two arbitrary constants, rather than one, which encode non-perturbative effects that are 
competing with the weak-coupling calculation at the considered accuracy. We have used the 
lattice data of Ref. |35j] to extract the value of the unknown piece of the three-loop singlet 
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potential. Our analysis indicates the following value of dsy. 

ds,s = i.nTol X 10^ (41) 

where the central value corresponds to the best fit of the N^LL curve. For those values, 
an excellent agreement with lattice data is achieved in the region where the weak-coupling 
calculation is rehable. 
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APPENDIX A: RELATION BETWEEN cq AND 03,^ 

Writing the expansion of the singlet potential (for /j, — 1/r) as 

V. = ^^f:v}-^al\lM, (Al) 

^ n=0 

we have 

V}'^ = (co + 27EC1 + (^47! + y) C2 + (87I + 27r'^E + 16C(3)) , (A2) 

where ior Uf — 

1639 1331 , 

ci = 290.769, C2 = ^, C3 = -^. (A3) 

K3) 



Finally, the relation between d^^s and Vs is 

_a; 

'(47r)3 



Cp^ = T/(3)_ (A4) 
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